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Coupled Oscillator Systems Having Partial VT Symmetry 
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This paper examines chains of N coupled harmonic oscillators. In isolation, the jth os¬ 
cillator (1 < j < A^) has the natural frequency ujj and is described by the Hamiltonian 
^p‘j The oscillators are coupled adjacently with coupling constants that are purely 

imaginary; the coupling of the jth oscillator to the (j + l)st oscillator has the bilinear form 
ijXjXj-^i (7 real). The complex Hamiltonians for these systems exhibit partial VT symme¬ 
try; that is, they are invariant under i ^ —i (time reversal), Xj —Xj {j odd), and Xj Xj 
{j even). [They are also invariant under i —i, Xj Xj {j odd), and Xj —Xj {j even).] 
For all N the quantum energy levels of these systems are calculated exactly and it is shown 
that the ground-state energy is real. When ujj = 1 for all j, the full spectrum consists of a 
real energy spectrum embedded in a complex one; the eigenfunctions corresponding to real 
energy levels exhibit partial VT symmetry. However, if the ujj are allowed to vary away 
from unity, one can induce a phase transition at which all energies become real. For the 
special case N = 2, when the spectrum is real, the associated classical system has localized, 
almost-periodic orbits in phase space and the classical particle is confined in the complex- 
coordinate plane. However, when the spectrum of the quantum system is partially real, the 
corresponding classical system displays only open trajectories for which the classical particle 
spirals off to infinity. Similar behavior is observed when N > 2. 

PACS numbers: 11.30.Er, 03.65.Db, ll.lO.Ef, 03.65.Ge 


I. INTRODUCTION 

There are many experimental and theoretical studies of T^T-symmetric coupled-oscillator Hamil¬ 
tonians [IHS]- In most cases the starting point is either a coupled set of 'PT-symmetric equations 
of motion, or a PT-symmetric Hamiltonian that governs such a system. It has been established 
that such PT-symmetric systems exhibit a rich phase structure with phase boundaries depending 
on the number of oscillators, how they are coupled, and the values of the coupling parameters [5]. 

In a recent paper on radiative coupling and weak lasing of exciton-polariton condensates, Aleiner 
et al. [7] considered a Hamiltonian function that governs condensation centers that are bilinearly 
coupled by a term of the form igzz"^, where each center is described by the complex coordinate z 
and ^ is a coupling strength. They investigated the classical dynamics of the system. While there 
is no obvious underlying symmetry, the authors found closed paths in their spin trajectories. This 
intriguing result motives the current study of an unusual type of oscillator system; namely, a chain 
of N harmonic oscillators with pure imaginary coupling. The Hamiltonian for the jth oscillator 
(1 < j < N) has the form ^p‘j + where the natural frequency uuj is real and positive. The 

jth oscillator is coupled to the {j + l)st oscillator by an imaginary coupling constant zy, where 7 
is real and independent of j. The coupling term is bilinear] that is, it has the form i^XjXjj^i. The 
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Hamiltonian that governs this system of N adjacently coupled oscillators has the form 


^ N N-l 

= 9 (d + m (^ > 2). 


J = 1 


J = 1 


( 1 ) 


This complex Hamiltonian is not VT symmetric because i changes sign under time reversal T 
and it is assumed that every coordinate xj changes sign under parity V. However, is partially 
VT symmetric; that is, it remains invariant if we change the sign of i and simultaneously reverse 
the sign of only the odd-numbered or only the even-numbered coordinates. To illustrate, we define 
Vj as the operator that reverses the sign of Xj but does not affect any other coordinate. Then, 
H2 is partially VT symmetric with respect to ViT and also with respect to V2T. Similarly, is 
partially VT symmetric with respect to V1V3T and also with respect to V2T. Note that reversing 
the signs of an even number of coordinates is achievable by a rotation but reversing the signs of 
an odd number of coordinates is not achievable by a rotation. For example, for = 2, ^ —xi, 

X 2 —X 2 is merely a rotation by an angle of tt in the xi^X 2 plane, but xi —xi, X 2 X 2 cannot 
be achieved by a rotation. For = 3, ViVs is a rotation but V2 and also V1V2V3 are not. 

Systems having partial VT symmetry have remarkable properties. In Sec. we set Uj — 1 for 
all j and show that for small N and for all values of the coupling parameter 7 the ground-state 
energy of the quantum system is real and positive. Then, in Sec. in we present the exact solution 
for the complete quantum spectrum for all N. We find that the ground-state energy is always real, 
but that the full spectrum is partly real and partly complex. For each energy, we calculate the 
corresponding eigenfunction and demonstrate that simultaneous eigenfunctions of the Hamiltonian 
and the partial VT operator have real energies, while those that are not partially VT symmetric 
are associated with complex energies. Thus, partial VT symmetry is associated with a partially 


real energy spectrum. In Sec. IV we relax the constraint that ujj = 1. We show that for V = 2 it is 


possible to choose the natural oscillator frequencies to make the energy spectrum completely real. 
Thus, there is a phase transition from a partially real to a completely real spectrum. This result is 
shown to hold in a modified form for A^ = 3 and N = 4. Next, in Sec.|V} we investigate the classical 
solutions for the N = 2 and N = 3 systems and find no remnant of the partially 'PT-symmetric 
phase; that is, all classical orbits are open unless the quantum spectrum is entirely real, in which 
case the orbits are all closed and periodic. Brief concluding remarks are given in Sec. |VI[ 


II. GROUND-STATE ENERGIES OF N COUPLED OSCILLATORS WITH ujj = 1 

In this section we show that the ground-state energy of a quantum system of N coupled oscillators 
with natural frequency ojj = 1 is real and positive. 


A. Two coupled oscillators 

Let us consider the quantum-mechanical Hamiltonian of two coupled oscillators H 2 — \p^ + 
-h \x‘^ -h \y‘^ + i"fxy^ where x and y are the coordinates, p and q are the conjugate momenta, 
7 is a coupling strength, and uoi = UJ 2 — This Hamiltonian is partially VT symmetric because 
it is invariant under the transformations VxT or PyT, where Vx : {x^y) (—x,^), Vy : {x^y) 

(x, —^), and T : i ^ —i. The Schrodinger equation associated with H 2 is 

- Idy + + i-fxy) 'tp{x,y) = Eip{x,y). (2) 

The ground-state eigenfunction has the (non-nodal) gaussian form 

'ipo{x^ y) — exp (—— \ay^ + bxy) , 


( 3 ) 
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where a and b are constants. Note that 'ipoix, y) is VT symmetric in either x or y. Inserting @ into 
Q and matching powers of x and y gives the three equations Eq — and 2ab = —i'j. 

The physically acceptable solution to these equations requires that b be imaginary, b = 
and that Eq = a he the real and positive solution to — 7^/4 = 0, 

Eq = a = + + 7^) • (4) 

Note that because b is imaginary and a is real and positive, ^o{x^y) vanishes as x‘^ + y‘^ ^ oc. 


B. Three coupled oscillators 

For three oscillators the Hamiltonian Hs in Q with uj = 1 has the form 

Hs = + ^x‘^ + ^y‘^ + + h{xy + yz). ( 5 ) 

Again, is partially VT symmetric; it is invariant under VyT (and also VxVzT). The lowest- 
energy eigenstate has the form ' 0 o(^? 2 /? z) = exp [—+ z^^) — ^by‘^ + c{xy + yz) + dxz'\ , where 
a, 6 , c, and d are constants. Solving the Schrodinger equation Hsipo{x,y^ z) = E^o{x^y, z) and 
comparing powers in x, and z gives the five equations Eq = a+^b^ 1 = + c^, 1 = 2c^ + 6 ^, 

iy = c{d—a—b), 2ad = c^. We solve these equations and verify that the eigenfunction is normalizable 
[^o{x^y^z) vanishes as x‘^ + y^^ + z^^ ^ oc] and that, even though is complex, the ground-state 
energy is real and positive, 

Eq — \ ~h \ {2 2^Jl + 27 ^^ . ( 6 ) 

The ground-state eigenfunction has partial 'PT-symmetry. Also, in the limit 7 ^ 0 the 

oscillators decouple and we recover the expected result that Pq = 3/2. 


C. Four coupled oscillators 


For four coupled oscillators the coordinates are x^y^z^w^ the canonical momenta are 
the Hamiltonian H 4 with cjj = 1 is partially VT symmetric in the variables x^ z or y, rc, and reads 

^4 = ^ {p^ + + 5^) -h ^ Ey‘^ E z‘^ E w‘^) -h i'j{xy + yz-h zw). (7) 


We solve the Schrodinger equation 7^4'0o = £’'00 with the ansatz for a partially PT-symmetric 
ground-state wave function of gaussian form 


^o{x,y,z,w) = exp 





-h c{xy + zw) -h d{xz -h yw) -\- exw -h fyz , 


( 8 ) 


where a, 6, c, d, e, / are six arbitrary constants. This leads to the conditions Eq = a + b^ p + d? + 


-\- b‘^ = 1, -h d^ -h 


-h = 1, iy = 2cd — 26/, cf + ce = bd -\- ad^ i'y = df -\- de — cb — ca, 
ae = cd. Clearly, the complexity of the coupled nonlinear system of equations increases rapidly as 
the number of coupled oscillators increases. 

For the case of N coupled oscillators with Uj — we show in Sec. HI that the ground-state 
energy Eq is 

N 


s/l + 2^7 COS [j7T/{N + 1)]. 


J=1 


(9) 
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By setting N = 2 or N = 3 , we readily recover Q and Q. For = 4 , Q yields the value 


En = 


/-1 r /- 

5 + 5 v ^ + ■v/ 5)/2 + ^ + 5^1 + 72(3 - a/ 5)/2 


1/2 


( 10 ) 


Closer inspection of Q reveals that the ground-state energy of such coupled oscillators is always 
real. Indeed, @ can be rewritten as 

^0 = ^ 5A/l-0 472cos2[j7r/(Al-0l)]) ^ . (11) 

J = 1 


III. EXACT EIGENFUNCTIONS AND SPECTRA OF N COUPLED OSCILLATORS 


A. Two coupled oscillators 


Let us return to the two-coupled-oscillator system governed by the Hamiltonian H 2 with cjj = 1. 
The transformation xi = {x + y)l\/ 2 , X2 — {x — y)/V^ decouples the oscillators, leading to the 
Hamiltonian H — ^pl -h ^(1 + i^)x\ + ^P 2 + ^(1 — which has complex-conjugate frequencies 

= 1 -h iy and ~ '^y* Apart from a normalization constant, the eigenfunctions are 

'^ni,n2{xi,X2) ^ (12) 

with corresponding energy eigenvalues £’ni,n2 = {pi + \) +'^2 (^2 + ^)' terms of the coupling 
parameter y the frequencies are 

i/i = I/* = (i lyrvyj ^ -f i -1- ^, 

whose real parts are positive. The general result for the energy spectrum is 

^ 711,712 ~ + y^^ (^1 + 712 + 1) + i + ^\/l + 7 ^^ (^1 ~ '^2)- 


Note that the spectrum is real if ni —1x2. If ni — 1x2 — 0 , we recover the ground-state energy 
in In addition, we obtain the corresponding eigenfunction from ( [T^ , 


'^Qfl{x,y) = exp 


\+ 




r 


+ y^)-i[-. 


+ yi + -iA 


1/2 


xy 


which verifies the ansatz © and explicitly demonstrates that an eigenfunction having the partial 
VT symmetry of the Hamiltonian is associated with a real eigenvalue. Note also that the real 
spectrum is part of a larger spectrum containing complex-conjugate pairs. This can be illustrated 
by the choice ni = 1 and 77.2 = 0 or ni = 0 and 77.2 = 1 : 


^1,0 


(2+2^1+72) i (^-1 +1 yr+y) 


and y) = \/ 2 (l + i'y)^^^{x -\- which is neither VxT nor 'PyT-symmetric. In addition, 

Po,i = ^10 ^o,i(^? 2 /) = ^i,o(^5 2 /)*- The real parts of the energies are 77.1 -h 77.2 + 1 -fold 

degenerate, as shown in Fig. 

The nature of the eigenfunctions associated with the first few energy levels is depicted in Fig. 
The ground-state (77.1 = 77.2 = 0 ) and the third-exited-state (77.1 = 77.2 = 1 ) eigenfunctions are par¬ 
tially VT symmetric, as can be seen in the left upper and lower diagrams, while the eigenfunctions 
corresponding to the first and third (complex) eigenvalues (77.1 = 0 , 77.2 = 1 and 77.1 = 0 , 77.2 = 2 ), 
which are not partially VT symmetric are shown on the right diagrams. 
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Rje(E) 



Im(E) 


State 


-4- State 


10 


Figure 1 : Real parts (left panel) and imaginary parts (right panel) of the first eleven energy levels of H2 
(TV = 2 , 7 = 1 ). 



Figure 2 : (Color online) Absolute value of the eigenfunction for the ground state (upper left), ni = 722 = 1 
(lower left), ni = 0 , 77.2 = 1 (upper right) and 77.1 = 0 , 77.2 = 2 (lower right) for 7 = 1 . 


B. Three coupled oscillators 

To find the exact solution to the Schrodinger equation for in ^ with Uj = 1 we make the 
transformation xi = {x — z)l v^, X2 = vlV^ + {x + z)/2^ X3 = —y/\f2 + (x + 2;)/2, which decouples 
the oscillators, giving H = + ^P2 + ^^2^2 + \P3 + with = 1 , z/l = 1 + 

and z/3 = 1 — Thus, the unnormalized eigenfunctions are 

'^{xi,X2,X3) = Mniixi)y.n2{V^X2)y-n3{V^X3) exp [-{xj + l' 2 xl + 










6 


and the energies are E = ni + ^ + 1^2 (^2 + ^) + ^3 (^3 + ^) 5 where 

^2,3 — + ^\/l + 27^^ ± i ^ + ^\/l + 27^^ 

Thus, the energy spectrum can be expressed as 

E = Til + ^ + \ l \ E \ \/^ E 2 ^ + \ l \ E \ '\/^ E 2 ^ (77.2 + 77.3) + — \ + ^ V^l + 27^ (77.2 — 77.3). 

Evidently, if the second and third oscillators are in the same state (77.2 = 77.3), the energy is real 
and the corresponding eigenfunctions are partially VT symmetric. In particular, the ground-state 
energy (§ is recovered and the ground-state eigenfunction is 


= exp -(1 -h a){x^ -h z^)/A - ay^/2 - (a - l)xz/2 - ib{xy -\- yz)lV2 
where a = Rez /25 b = Imz/ 2 - The first ten energies of for 7 = 1 are shown in Fig. 


Rje(E) 

4 - 


Im(E) 


State 


State 


Figure 3 : Real parts (left panel) and imaginary parts (right panel) of the first ten energies for (7 = 1, 
N = 3). 


C. Four coupled oscillators 

For the Hamiltonian 0 , which governs four linearly coupled oscillators, the transformation 


Xl,2 = 


a^ 3,4 = 


1 


2V5 . 
1 


2^5 


\Jb — Vb (x ± ru) ± \Jb + Vb {y ± z) 
\Jf> + \/5 {x^w):k \J— \/5 (y ^ z) 


exactly decouples the oscillators. The new Hamiltonian takes the form 

TT 12 | 122 | 12 | 122 | 12 | 122 | 12|122 

H^2Pi + 2^1^1 + 2P2 + 2^2X2 + 2^3 + 2^33^3 + 2Pa + 2^aXa^ 

where — U 2 * = 1 + 5 ^ 7(1 + \/5) and r'f = r'l* = 1 + ^i 7 (—1 + \/5) are complex frequencies. Let 
Re ui = Re z/2 = ^5 Re ^3 = Re 7/4 = C, Im ui = —Im z/2 = Im z/3 = —Im z/4 = T), where 

\j±\ + \^i + 57^(3 + a/ 5), (c,r>) = 
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In terms of these variables and the quantum numbers ni, n 2 , ns, n 4 , the total energy is 


^711,712,713,714 — ^{'^1 + n2 + l) + C + n4 + l) + iB[ni — 112 ) + iD{n^ — n^). 


Thus, if ni = n 2 and ns = n 4 , the energy is real. When the energy is real, the correspond¬ 
ing eigenfunction is always VxzT- or "PywT- symmetric. For example, the ground-state energy is 
£’o,o,o,o — A-\- C and the corresponding eigenfunction is 


^o,o,o,o = exp 


--(^+c+ 


~{A - C){xz + yw) - ^ 



i fr. B ^ D\ 


i 

7 ! 


{B -h D){xy -h wz) 


(13) 


This eigenfunction displays the symmetries assumed in the ansatz Q. 

As another illustration, we consider the case in which the first two oscillators are in the first 
excited state, and the other two in the ground state (ni = 77.2 = 1, ^3 = 77.4 = 0). The energy is 
E’l, 1 , 0,0 = 3A -h C and the eigenfunction is 


^ 1 , 1 , 0,0 = -y^A^^B‘^ (5 — — uP) — (5 -h V^){y^ — z^) + 4^/5xz — A\/byw 

o 

Once again, the energy is real and the eigenfunction is VxzT- or VywT-symmeiAc. 
energy 0 , 0,0 = 2A + C + iB arises for the choice 77.1 = 1, 77.2 = 77.3 = 77.4 = 0. 


^o,o,o,o- 
A complex 


D. Five coupled oscillators 


The Hamiltonian for five coupled oscillators is 

= \ + 5^ -h ^ -\-y‘^ -\- z‘^ + w‘^ -\- u^) + iy{xy yz + zw + wu). 

Rather than decoupling the oscillators we treat this case by constructing the secular equation 

det [Mjk - u^Sjk) = 0, (14) 

where Mj]^ is the tridiagonal matrix defined as 


Mjk — 


d^U 


dqjdqk 


/ 1 ij 0 0 0 \ 

iy 1 iy 0 0 

0 iy 1 iy 0 

0 0 iy 1 iy 

Vo 0 0 iy 1 y 


(15) 


U is the potential, and qj and Qk are coordinates. 

The solution to the secular equation (14) gives complex-conjugate pairs of frequencies and one 
real frequency: 2 = 1 ± ^yv^, ^ 3^4 = 1 ± iy, = 1. Thus, the decoupled Hamiltonian is 

Tj 12|122|12|122|12|122|12|122|12|12 
= 2^1 + 2 ^ 1^1 + 2P2 + 2 ^ 2^2 + 2^3 + 2 ^ 3^3 + 2^4 + 2 ^ 4^4 + 2^5 + 2^5 

and the energy of the system reads 

E = ui (ni + ^) + z/2 ( 77.2 + ^) + ^3 (^3 + ^) + ^4 (^4 + ^) + ^5 + ^5 


where 77 . 1 , 77 . 2 , 77 . 3 , 77 . 4 , 77.5 are nonnegative integers. 














E. General case: N coupled oscillators with ujj = 1 


In this section we consider the Hamiltonian Q for N linearly coupled oscillators with ojj = 1 . 
To obtain the frequencies of the decoupled oscillators we use ( 14 ) to construct the N xN tridiagonal 
matrix secular equation, det(M — z/^I) = 0 , which has the form 


Dn = 


1 — iy 

l-z/2 


Z7 


I'j 


iy 1 — 


i'y 1 — iy 
i'j 1 — z/^ 


= 0 . 


Because this matrix equation is tridiagonal, satisfies the three-term recurrence relation 
Dk + - l) Dk-i - l‘^Dk-2 = 0 (k = 1 , 2 ,... ,N), 

where = 1 and D-i — 0 . We solve this difference equation to obtain the frequencies 

z/^ = 1 + 2 i'j cos[j7T/{N + 1 )] (j = 1 , 2 ,..., N). 


Thus, the exact expression for the total energy of the system of N oscillators is given by 

N 

E = yi + 2i7COs[j7r/(iV + 1 )] {rij + i) , 

where nj > 0 {j = 1 ,..., A^). Choosing nj = 0 for all j, we find the exact ground-state energy in 
( [TT| ), which has been shown to be real. 


IV. COUPLED OSCILLATORS WITH ARBITRARY FREQUENCIES 

In Secs. and |lll| the oscillator frequencies multiplying xj were set to unity. Our conclusion 
in the foregoing analysis was that a real spectrum is embedded in a complex spectrum containing 
complex-conjugate pairs of energies. We now relax this constraint on the natural frequencies. For 
the two-, three-, and four-coupled-oscillator systems, we demonstrate that for an appropriate choice 
of (jjj the spectrum can be entirely real. 


A. Two coupled oscillators with general natural frequencies and ujy 

The Hamiltonian H2 in (|l| reads H2 = + \(f + + \^yy‘^ + i"fxy. The frequencies of 

the decoupled oscillators in this case are 

^1,2 = + ^{ul - ujIY - 472j ( 16 ) 

and the energies of the system are (^1 + ^) + ^2 (^2 + where ni, n2 > 0. 

In contrast to the results found in Sec. m the entire energy spectrum can be real for specific 
values of ojx^ and 7. For this to be so, the parameters must satisfy the condition 

> 2 | 7 |. 


( 17 ) 
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The case considered in Sec. Ill had ujx = ooy = I and 7 = 1 , which does not satisfy this condition 


and, as we saw, the energy spectrum was only partially real. 

To have real eigenvalues the associated eigenfunctions must all be partially VT symmetric. This 
can be seen explicitly by decoupling the oscillators with the transformation 


xi = y/A-\- B{Dx + Ex + iCy)/{ 2 CE ), X2 = VA — B{—Dx + Ex — iCy) / { 2 CE ), 

where A = 87^ — — Wy)‘^, B = — Wy)[{w‘^ — WyY — 47^]^^^, C — 27, D = w: 

E — — WyY — 47^] leading to the Hamiltonian H = ^pl + ^i^ixf + ^P2 + ^ 

and z /2 are given in ([T^. Up to a normalization constant, the eigenfunctions are 


'^ni,n2{xi,X2) = 'Hni {\/i^Xi)'Hn2 (^1^X2) ex.p - ^^2x1) . 

Note that A -\- B and A — B have opposite signs in the T^T-symmetric phase and that A — B < 
0 . Rewriting 4^711,712(^15^2) in terms of the original variables x and y, one can show that the 
eigenfunction has partial VT symmetry because 

(7^x7^)4^711,712 (^5 2/) ~ (~1) ^ ^771,772 (^5 2/)? (7^^7^)4^771,772 (^5?/) ~ (~l) ^4^771,772 (^5?/)- 


To illustrate, we consider the case — 3 , ojy = 1 , and 7 = 1 / 2 . The relation ( 17 ) is satisfied 
and we find the purely real nondegenerate spectrum shown in Fig. When — cjj < 2|7|, the 
energy spectrum is only partially real. Thus, there is a phase transition from the unbroken partially 
PT-symmetric phase to the broken one. For example, keeping Uy — 1 and 7 = 1 / 2 , but adjusting 
uj‘^ so that it passes 2, the first-excited-state energy becomes complex. 

7 -| 

E 

6 - _ 

5 - “ 

4 - “ 

3 - — 

2 - 
1 - 
0 - 


Figure 4 : First ten energies of H2 for the parameter choice = 3 , = 1 , and 7 = 1 / 2 . These states have 

the quantum numbers (711,712) = ( 0 , 0 ), ( 0 , 1 ), ( 1 , 0 ), ( 0 , 2 ), ( 1 , 1 ), ( 0 , 3 ), ( 2 , 0 ), ( 1 , 2 ), ( 0 , 4 ), ( 2 , 1 ). 


B. Three coupled oscillators with general natural frequencies ujy and ujz 

For the three-oscillator Hamiltonian + + + ^(jjyy‘^ + i'j{xyEyz) 

the frequencies = A of the decoupled oscillators satisfy the cubic equation /(A) = 0 , where 

/(A) = V - + ul + ul) A 2 + {ulul + + 27^) A - ^yy - y + 7 ^ 

If the discriminant associated with this equation is positive, three real distinct roots can emerge, 
giving a real spectrum. To guarantee that the roots are positive, it is necessary that 

/(O) < 0, Amax ^ 0, Ajnin ^ 0, /(A max) ^ O5 /(Amin) < 0 
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are all fulfilled, with Amin and Amax being the extrema of the polynomial. Figure displays the 
regions in which the frequencies of the decoupled oscillators are all real (blue shaded areas) in the 
parametric space of Ux and Uz for fixed values of 7 and ujy. This figure shows that several regions 
of unbroken symmetry exist. This is in contrast to the case of the two coupled oscillators. 



Figure 5 : (Color online) Unbroken partial T^T-symmetric phases of depicted as shaded areas for the 
parameters 7 = 1/12 and = 2/3 (left panel) and 7 = 1/3 and ujy = l (right panel). 

For example, Figure shows that uj‘^ = 1 / 3 , Uy — 2 / 3 , = 1 , and 7 = 1/12 gives an unbroken 

symmetry phase. We obtain three different real positive (decoupled) frequencies 

I/I = 773, z/2 = \/(8 +Vl4)/12, z/3 = 

Thus, the spectrum is entirely real with energies given by Eni^n2,nz — (^1 + ^) + ^2 (^2 + ^) + 

^3 (^3 T ^)- By fixing only 7 we can find regions in the three-dimensional parameter space of ujx^ 
(jjy^ and (jJz for which unbroken symmetry (and therefore a real spectrum) exists. This is shown in 
the colored volumes depicted in Fig. |^for the specific choice 7 = 1 / 12 . 


C. Four coupled oscillators with general natural frequencies ^nd 


The previous analysis can be applied to the four-coupled-oscillator Hamiltonian 

TT 12 | 12 | 12 | 12 | 122 | 122 | 122|12 2 / , , \ 

-h -\-^ujx^ ^ Ei^[xy Eyz + zw), 


where cjy, and Uyj are real frequencies and 7 is a real coupling parameter. 
The eigenvalues of the matrix 


/ n 0 0 \ 

i7 i'y 0 

0 i7 7 *7 

\ 0 0 ij ojI J 


M = 
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1.0 



1.0 


Figure 6: (Color online) Unbroken partial T^T-symuietric phases of for 7 = 1/12 depicted as colored 
voluuies. 


are the squares of the corresponding decoupled-oscillator frequencies. The eigenvalues — X satisfy 
the fourth-order equation /(A) = — aX^ + bX^ — cA + d = 0 , where 




''X ' 

L 22, 22, 22, 22, 22, 22io2 


y'^w 


222, 222, 222, 222io22io22i 22, 22 

C = + ^x^z^w + ^y^z^w + 27 + 27 + 7 + 7 , 


''x'^y'^w 


j 2222 , 222 , 222 , 222 , 4 

^ ^x^y^z^w T ^x^y T ^x^w ^z^w T * 


Regions in which all decoupled oscillator frequencies are real give a completely real energy 
spectrum, which means that partial VT symmetry is unbroken. This requires that /(A) have four 
positive roots, which is the case if /(O) > 0 . In addition, if f'{X) has three positive roots, the 
extrema of /(A) lie on the positive abscissa. To have four real roots the minimum value of /(A) 
must be negative, and the maximum value must be positive. Figure shows the regions in which 
these conditions are fulfilled (that is, the regions in which the partial VT symmetry is unbroken) 
for specific choices of the parameters. Fixing the values of = 1 and = 4 as in the right panel 
of Fig. we can investigate the development of the phase boundaries as a function of the coupling 
strength 7. This is shown in Fig. 


V. CORRESPONDING CLASSICAL THEORY 

In this section we investigate the features of partially T^T-symmetric classical theories. We begin 
with the two-coupled-oscillator Hamiltonian H2 with frequencies = cjy = 1 . Hamilton’s classical 
equations of motion lead to 

x\t) + x{t) + ijy{t) = 0, y\t) + y{t) + ijx{t) = 0 
and combining these equations gives the fourth-order differential equation 

x""{t) + 2 x"{t) + (1 + 7^)x(t) = 0. 
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Figure 8: (Color online) Unbroken partial PT-symmetric phases of P4 in the (ujx^ooy) plane depicted as 
shaded areas for the parameter choices = 1 and = 4 for twenty values of 7. 
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We seek solutions x{t) = and find that A = satisfies the quadratic equation A^ — 2 A + l+7^ = 
0 , so z/ = ±^1 zb i'j. Thus, the characteristic frequencies are always complex. By decomposing 
y/1 + iy = a -\- ib into its real and imaginary parts we can write the general solution as 


x{t) — {A-\-D)e + {B + C)e^^ cos(at) + i {A — D)e {B — C)e^^ sin(at), 








where A, S, C and D are arbitrary constants. Therefore, for any initial conditions, the real 
and imaginary parts of x{t) are oscillatory and growing (or decreasing). As a consequence, the 
trajectories in the complex-x plane spiral outward (or inward). Hence, the classical paths are open 
(see Fig. left panel). Thus, although the quantum spectrum is partially real, this partial reality 
does not give rise to closed classical trajectories. 



Im(x) 



Figure 9 : Left panel: classical trajectory in the complex-x plane for H2 with uj^ = ujy = 1 with the initial 
conditions x{ 0 ) = y{ 0 ) = —1 — i and i:( 0 ) = ^( 0 ) = 1 + ^/2, and 7 = 1. Right panel: Classical trajectory in 
the complex-x plane for H2 with = 2, = 1, and 7 = 1/2 for the initial conditions x(0) = ^(0) = 

and x(0) = ^(0) = 2 — i/A. 

More generally, if the coupled oscillators described by H2 have natural frequencies uox and uoy^ 
we obtain the equation 


x"\t) + {w^ + u;^) x\t) + {w‘^Wy + 7^) x{t) = 0. 

Again seeking solutions x{t) = we find that z/^ = + ^Wy ± ^ [(cj^ — cj^)^ — We 

deduce that four real values of z/ exist when \uj‘^ — > 2|7|, which is precisely the condition that 

guarantees a fully real spectrum in the quantum system [see Thus, the transition from the 

broken partial-'PT-symmetric phase to the unbroken phase occurs at the same point as for the 
quantum case. For the parameter choice — 2 , ujy — 1 , and 7 = 1/2 the classical trajectory 
depicted in Fig. (right panel) spirals outward, which indicates that VT symmetry is broken even 
though the system is partially VT symmetric. The behavior of the trajectories in the unbroken 
phase is illustrated in Fig. Here, one sees that while the classical trajectory is not closed it is 
confined to a compact region in the complex-x plane. We thus observe the phase transition at the 
classical level; that is, we observe the transition from spirals (broken phase) to localized trajectories 
in the complex-x plane (unbroken phase), which happens at \u‘^ — Uy\ — 2\^\. 
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Im(x) 



Figure 10: Classical trajectory in the couiplex-x plane for H 2 with the parameter choice = 3, = 1, 

7 = 1/2, with the initial conditions x(0) = —1 + ^(0) = 2 — 2i, x(0) = 1 + i/2, ^(0) = 3/2 + i (left panel) 

and x(0) = 1 + 2i, ^(0) = 2 + i/4, x(0) = —3 + i, ^(0) = 1/2 + 5i (right panel). 


Rc(q) Rc(q) 



Figure 11: Poincare section of the classical trajectory in phase space for H 2 with the parameter choice = 3, 
= 1, 7 = 1/2 and initial conditions x(0) = 1 + i/2, ^(0) = —2 + 2i, i(0) = 1 + i/2, ^(0) = —3/2 — 3i/2 
(left panel) and x(0) = 2 + 5i/2, ^(0) = 2 — 3i, x(0) = 1 — i/2, ^(0) = 2 — 2i (right panel). The structure of 
this plot indicates that the orbits are almost periodic. 


In addition to studying the trajectory in the complex-x plane, one can study the trajectories in 
phase space by plotting a Poincare section. From the structure of the Poincare plot we conclude 
that the confined trajectories are almost periodic [8]. This is illustrated in Fig. 11 The appearance 
of open almost-periodic trajectories occurs because the number of degrees of freedom exceeds one; 
if = 1, the classical orbits associated with real quantum energies are closed [9]. 

A similar analysis can be done for the three-coupled-oscillator Hamiltonian with general 
natural frequencies. The classical equations of motion are 


x''{t) + ujlx{t) + i 7 y(t) = 0, y''{t) + + ii{x{t) + z{t)) = 0, z''{t) + ulz{t) + i'^yit) = 0. 

Seeking solutions of the form x{t) = y{t) = and z{t) = we obtain a cubic 
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Im(x) 



Im(x) 



Figure 12: Classical trajectory in the complex-x plane for Hs with the parameter choice = 1/3, ujy = 
2/3,7 = 1/12, = 1 (left panel) and = 13/20 (right panel) with the initial conditions x(0) = —2 + i, 

^(0) = 3 — 3i, 2:(0) = 3 + 2i, x(0) = —1 + 3i, ^(0) = 3 + 2i, i:(0) = —2 + i. 


Rc(q) 



Figure 13: Poincare section of the classical trajectory for Hs with the parameter choice = 1/3, = 2/3, 

= 1, 7 = 1/12 and the initial conditions x(0) = —2 + i, ^(0) = 3 — 3i, ^(0) = 3 + 2i, j:(0) = —1 + 3i, 
^(0) = 3 + 2i, i:(0) = —2 + i. This figure indicates that the orbit is almost periodic. 


equation for A 




A 3 - (c.2 ,,2) ^2 + (^2^2 + ^2^2 + ^2^2 + 2^) A - ^ 7^ = 0. 

The characteristic frequencies are real if and only if the corresponding quantum system is in an 
unbroken PT-symmetric phase (all eigenvalues are real); the criteria for real eigenvalues is given 
in Sec. 0 To illustrate, recall that in Sec. 


IV 


the parameter choice = 1/3, ojy = 2/3, = 1, 

7 = 1/12 lies in the unbroken phase and gives a real energy spectrum. Figure shows that for this 
parameter choice, the classical trajectory is confined to a compact region in the complex-x plane, 
whereas for the choice = 1/3, Uy = 2/3, = 13/20, and 7 = 1/12, for which the quantum 

symmetry is broken, the classical trajectory spirals outward to infinity. A Poincare section is given 


in Fig. 13 for the parameter choice of Fig. left panel. 
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VI. BRIEF CONCLUDING REMARKS 

In this paper we have examined systems of N linearly coupled oscillators that are partially VT- 
symmetric. In the quantum-mechanical analysis we have found that the ground state of each of 
these systems is always real. We have shown that the entire spectrum may in fact be completely 
real depending on the values of the natural frequencies ojx^ • • • and their relation to the 

coupling strength 7 . This happens even though the system is only partially VT symmetric. We 
have studied this in detail for systems of two and three coupled oscillators. A phase transition point 
exists beyond which the energy spectrum is only partially real. 

For the two and three classical oscillator systems, we find a phase transition at exactly the same 
point as the quantum-mechanical oscillator systems. When the eigenvalues of the quantum system 
are all real, the classical trajectories are confined and almost periodic, but when the quantum 
eigenvalues are partly real and partly complex, the corresponding classical system always has open 
trajectories that spiral out to infinity. 
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